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Abstract 

We present a new approach for studying the problem of optimal 
hedging of a European option in a finite and complete discrete-time 
market model. We consider partial hedging strategies that maximize 
the success probability or minimize the expected shortfall under a cost 
constraint and show that these problems can be treated as so called 
knapsack problems, which are a widely researched subject in linear 
programming. This observation gives us better understanding of the 
problem of optimal hedging in discrete time. 

1 Introduction 

In this paper we work with a finite, arbitrage-free and complete market 
model in discrete time. It is a well-known fact that a European option 
can be hedged perfectly in such a market and that the price of a 
perfect hedge is equal to the unique arbitrage-free price of the option. 
However, the seller of the option might not be willing to use all the 
money he has received from selling the option to construct a hedging 
strategy. She/he may instead want to create only a partial hedge. This 
kind of a hedge costs less than a perfect one but forces the investor to 
face the risk of a shortfall. 

There are different ways to measure the risk that an investor who 
uses a partial hedge must take. In this paper the optimality is mea- 
sured in terms of success probability (i.e. the probability that short- 
fall will not occur) and in terms of expected shortfall, when a cost 
constraint is given. Follmer and Leukert [4], [5] and Follmer and 
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Schied [6] have studied these problems in both discrete and continuous 
time. They search an optimal solution among admissible strategies, 
i.e. strategies which are self-financing and whose value processes are 
non-negative. Their solution techniques are mainly based on different 
applications of Neyman-Pearson lemma. Runggaldier, Trivellato and 
Vargiolu [10], Favero [2], Favero and Vargiolu [3] and Scagnellato and 
Vargiolu study the problems via dynamic programming in bino- 
mial and multinomial models in the case where the strategies are only 
required to be self-financing. 

The main contribution of this paper is that these types of partial 
hedging problems can be reduced to knapsack problems (see Section 
12.31 for a description of a knapsack problem). In particular, this new 
approach allows us to prove some of the existing results in an alterna- 
tive way. Gathering the results under knapsack theory helps us also 
to see interestingly how admissibility condition affects the optimal so- 
lutions. 

The partial hedging problems covered in this paper are presented 
in Section [2l In Section [Sj it is shown that the problem of maximizing 
success probability under admissibility constraint can be reduced to 
a 0-1 knapsack problem. An approximative algorithm for solving the 
problem is obtained by studying a related continuous knapsack prob- 
lem. Comparison with the Neyman-Pearson-based results in Follmer 
and Schied ISj is also carried out. In Section |4] we will prove that if ad- 
missibility constraint is omitted, the success probability is maximized 
by an optimal quasi-replicating strategy, i.e. a strategy that replicates 
the option in all but one state. In Sections [5] and [6] we consider the 
problem of finding a strategy that minimizes the expected shortfall. 
A knapsack problem approach is used to show how admissibility con- 
straint affects the solutions. An alternative proof for the result in 
Scagnellato and Vargiolu [H], p. 148, will be given. 

2 Partial hedging problems and knapsack problems 

2.1 Market model 

The market model that we work with in this paper is based on a finite 
filtered probability space {Ti}J^Q,F), where i7 = {ui, . . . ,oJn}, 

= {0,^7}, Tt = T = V{n) and F{{LOi}) > for all iOi G Q. Here 
V{0,) denotes the set of all subsets of Q. We set T to be equal to the 
maturity date of the European option which we want to hedge. 
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The prices of the d + 1 assets on the market follow a d + 1-dimen- 
sional non-negative {^j}-adapted stochastic process {St,t = 0, . . . , T}, 
where 

(i) 

and 5j is the value of asset i at time t. However, we will follow the 
approach in Follmer and Schied [6] and present all values in units of 
the numeraire asset 5*^, whose value therefore is assumed to be strictly 
positive at all times. The discounted price of asset i at time t (i.e. its 
value in units of the numeraire asset S^) is given by X^^^ := sj.^^ / sj:^^ 
and the corresponding price process is {Xt,t = 0, . . . , T}, where 

A trading strategy is an M'^+^-valued {^i}-predictable stochastic 
process ^ = {^^, t = 0, . . . , T}, where 

(i) 

and Q is the quantity of asset i in the portfolio at time t. The value 
process V = {Vt, i = 0, . . . , T} of ^ is defined through 

Vt=-C,-x, = Y,cPx^'\ 

i=0 

A strategy ^ is called self-financing if 

Xt = Ct+i ■ 

for alH = 0, . . . , T — 1. A self-financing strategy ^ is called admissible 
if its value process satisfies Vt > for alH = 0, . . . , T. 

We assume that the market is arbitrage-free and complete and let 
P* denote the unique equivalent martingale measure. Moreover, H is 
a discounted European contingent claim, i.e. a non-negative random 
variable on (r2,.FT,P). It is known that H can be hedged perfectly, 
i.e. there exists a self-financing strategy ^ such that H = S^rp- Xt- The 
initial cost of this replicating strategy is equal to the unique discounted 
arbitrage-free price of H and is given by vr^ = W{H), where vr^ is 
expressed in units of the numeraire asset. 

Remark 2.1. Completeness is usually defined to mean that every non- 
negative J^T-tne.asurahle random variable can he replicated. However, it 
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is easy to show that in our finite market the same holds for an arbitrary 
random variable. This can be done by studying separately its positive 
and negative part. 

Remark 2.2. It is well known that the value process V of a self- 
financing strategy ^ satisfies 

t 

Vt = Vo + Y,lr-{^i-^i-l)- 
i=l 

Moreover, in a finite probability space it is easy to show that V is 
always a ¥* -martingale. 

2.2 Partial hedging problems 

Next we fix an upper bound v < K*{H) for the initial payment (in 
numeraire units) that the investor wants to use for creating a hedging 
strategy. Under this cost constraint we try to find an optimal partial 
hedging strategy. The problems considered in this paper are stated 
below. In problems A and B optimality is measured in terms of suc- 
cess probability, in C and D through expected shortfall. Note that in 
problems A and C the optimal solution is searched among admissible 
strategies, whereas in B and D we only require that the strategies are 
self-financing. 

Problem A. Find an admissible strategy whose value process V max- 
imizes ¥(Vt > H) under the constraint Vq <v. 

Problem B. Find a self-financing strategy whose value process V max- 
imizes ¥{Vt > H) under the constraint Vq < v. 

Referring to the discussion in Follmer and Schied [6], p. 341 we state 
that the problem of minimizing the expected shortfall E[(i7 — Vt)"*"] 
can be simplified to the following problems: 

Problem C. Find an admissible strategy whose value process V max- 
imizes E(Vt) under constraints Vt < H and Vq <v. 

Problem D. Find a self-financing strategy whose value process V 
maximizes E(Vt) under constraints Vt < H and Vq < v. 
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2.3 The knapsack problem 

The main purpose of this paper is to show that the partial hedging 
problems A, C and D above can be reduced to knapsack problems. 
The knapsack problem is usually illustrated as follows (see e.g Dantzig 
[1], p. 273 or Martello and Toth ^Sj, p.l): A traveller has to fill a 
knapsack of a certain size c by selecting some of n objects having 
sizes Wi, i = 1, . . . ,n, respectively. The "comfort" or "gain" given by 
the objects is measured with numbers gi, i = l,...,n, respectively. 
The traveller wants to select objects that give her/him the maximal 
total "comfort" or "gain" under the constraint that the total size of 
the chosen objects will not exceed the knapsack size c. We model a 
possible decision by an n-dimensional binary vector x whose elements 
satisfy 

_ J 1 if object i is selected 
* 1^ otherwise. 

Mathematically, we have to find an n-dimensional binary vector x that 
maximizes 

n 
1=1 

among all binary vectors satisfying 

n 

wiXi < c. 

i=l 

This problem is commonly referred to as the 0-1 knapsack problem. If 
it is possible to choose any fraction of an object, i.e. if the decision 
vector can be of the form < Xj < 1, i = 1, . . . , n, we call the problem 
a continuous knapsack problem. The solution value of a continuous 
knapsack problem is clearly an upper bound for the solution value of 
the corresponding 0-1 knapsack problem. 

3 Maximizing the success probability under admissibil- 
ity condition (Problem A) 

We will first present briefly how Problem A is treated using Neyman- 
Pearson lemmas in Follmer and Schied [6], pp. 333-339. After that 
it will be shown that the same results can be accomplished through 
knapsack problem theory. Finally, we will discuss an approximative 
algorithm for solving Problem A. 
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3.1 Neyman-Pearson approach 

The following result can be found with its proof in Follmer and Schied 
[6], p. 335. Recall that v is the upper bound for the initial payment we 
want to use to set up a hedging strategy. 

Theorem 3.1. Assume that the set G maximizes the probability 
P(r) among all sets T £ Tt satisfying the constraint 

E*{H-lr)<v. (1) 
—A 

Then the replicating strategy ^ of the option H := H ■ IpA solves 

Problem A. Moreover, = {V^ > H}, where is the value process 
—A 

of the strategy ^ . 

In Follmer and Schied [6J the authors define the measure 

H 



E*{H) 



dF* 



and consider the generalized density dF/dQ that is received from the 
Lebesque decomposition. For our purposes it is enough to know that 
in our finite probability space this density takes the form 

, V- / P^/^l^ '^^'^i) ^ (0) 

dQ^ '^-'X +00 if Q(a;0 = O, 



where 
and 



E*{H) 

Once a := v/W^H) is defined, the level 



Pi := FiiJ,) (3) 

F*{u;,)H{u^,) 
i^i) = ^=^7777^ • (4) 



inf <{ c G M+IQ ( ^ > c-E*(F) ) < a [> (5) 



is introduced and it is shown, using Neyman-Pearson lemma, that if 

q{dF/dq > c* ■ E*{H)) = a, (6) 

then = {dF/dQ > c* ■ E*{H)} is an optimal set described in The- 
orem 13.11 
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To overcome the case when equation ([6]) is not satisfied, the in- 
dicator function IpA is replaced by a randomized test, i.e. an J-t- 
measurable function ip such that < < 1- After defining TZ as the set 
of all randomized tests, the authors consider the optimization problem 
of finding a randomized test ip* ^ TZ that maximizes the expectation 
E('i/;) among all '0 £ 7^ satisfying the constraint W[H ■ ip) < v. Such 
an optimal randomized test is by generalized Neyman-Pearson lemma 
given by 

'^'^^ = ^{dF/dQ>c*-E*{H)} + 7 ■ l{dP/dQ=c*-E*(/f)}5 (7) 

where 

_ a - Q{dF/dq > c* ■ E*{H)) 
^ ~ Q{dF/dQ = c* ■E*{H)) ' 

Note that Q{dF/dQ = c* ■ E*{H)) > in our finite probability space, 
i.e. 7 is well-defined. 

Finally it is shown that the replicating strategy ^ of the option 
H* := H ■ ip* actually maximizes the expectation of the so called 
success ratio tpy among all admissible strategies with Vq < v and that 
the success ratio tpv* of ^ coincides with -0*. 



3.2 Knapsack approach 

In our finite market model the problem of finding an optimal set 
described in Theorem 13.11 is. in fact, a 0-1 knapsack problem. To see 
this, note that for any F G Tt we have that 

Thus, using ([4)) we can write the constraint ([I]) in form 

Further, P(F) can be written as Y^ujierP^' where pi is as in ([3]). 

Since $7 consists of n elements, we see that finding an optimal 
success set F is equal to finding an optimal n-dimensional binary vector 
x^. The problem of finding F"^ can thus be written as the following 
0-1 knapsack problem: 
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Problem A'. Find an n-dimensional binary vector that maximizes 

n 

1=1 

under the constraint 

n 

qiXi < a. 

i=l 

We will below use the notation for the value of the optimal 
solution, i.e. 

n 

z^:=Y,Pr4=nyT>H), (8) 

i=l 

where the latter equality follows from the fact F"^ = > H} men- 
tioned in Theorem 13.11 

Many numerical algorithms have been developed to solve 0-1 knap- 
sack problems. A nice overview of some of these techniques can be 
found in Martello and Toth [8] and in Martello et al. [9]. However, 
exact solution algorithms can be difficult to implement in a large prob- 
ability space, where n is huge. A good approximative algorithm will 
be given in Section 13.31 

Remark 3.2. It is often assumed (see e.g. Martello and Toth J^, 
p. 14), that the variables pi, qi and a in a 0-1 knapsack problem are 
positive integers. This assumption is a corner stone even for some nu- 
merical algorithms. However, in Problem A' above we allow < < 1 
and < < 1. Algorithms that are based on the integer assump- 
tion, cannot naturally be used in this case. Note that our assumption 
V < 'E*(H), or equivalently a < 1, rules out the possibility of the trivial 
solution Xi = 1 for all i = 1,2, ... ,n. 

Next we will show that finding an optimal randomized test ^p* de- 
scribed in section 13.11 can alternatively be seen as a continuous knap- 
sack problem (cf. Section [2?3]l . The key point is that in our finite 
probability space there is a one-to-one correspondence between the set 
of all randomized tests and the set of all n-dimensional vectors x satis- 
fying < Xj < 1 for i = 1, . . . , n. This correspondence is simply given 
by equation 

^p{uJi) = Xi. (9) 
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Using ([3|) and dlj we see that the problem of finding an optimal V'* is 
equivalent to finding an n-dimensional vector x* that maximizes 



i=l 

under constraints 

n 

QiXi <a:= ^^^jjy 0<Xi<l, i = l,...,n. 

Assume now that the array (wi, . . . , Wri) is ordered so that the quo- 
tient Pi/qi is non-increasing. Here we use the convention pi/qi = +00, 
if qi = 0. Thus, the most preferable states, the items that give the 
highest probability compared to cost, are placed first in the array. 
Suppose now that we consecutively choose the items, starting from 
the one giving the best probability over cost quotient and continuing 
until we find the first item s that we no longer can afford to choose. 
In other words, we define the critical element 



(10) 



Due to assumption a < 1 we know that 1 < s < n. In Martello 
and Toth [8], pp. 16, 17 it is proved, by a simple contradiction, that an 
optimal solution to the continuous knapsack problem is given by 




x^ = < 



1, i = 1, . . . , s — 1 

^ = s (11) 
0, i = s + 1, . . . ,n. 



We denote the optimal solution value by 

n 



n 

z- := > piXi- 

i=l 



The value z* is an upper bound, the so called Dantzig 's bound for the 
optimal value z^ in jS]). If X]i=i 9« — then z^ = z* . Note the 
analogy with the Neyman-Pearson based result in jG]). 

Note that the difference between the optimal solutions ip^^ in ([7} 
and ip*{uji) = X*, i = 1, . . . ,n (cf. (fTTjl ) is that the Neyman-Pearson 
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approach tells us to search a critical level L := {dP/dQ = c* • W{H)} 
whereas the knapsack approach suggests looking for a critical element 
ujg. The definitions of c* in ([5|) and s in ([TO]) give that {ujg} C L. In 
fact, there is no restriction on the values of an optimal randomized 
test ijj on the critical set L except that the level condition 

E^(V'l{dP/dQ=c*.E*(//)}) = a - Q(dIP/dQ > c* ■ E*{H)) 

must be satisfied. This fact is mentioned also in Follmer and Leukert 
[5],p.l26. 

The particular form in ([2]) that the derivative dF/dQ takes in a 
finite probability space allows us to interpret the Neyman-Pearson re- 
sult in compliance with the knapsack approach, i.e. that we choose to 
hedge against the states having the best probability over cost quotient. 
This "cost-effectiveness" interpretation of the Neyman-Pearson lemma 
is mentioned already in Kadane [7j. 

3.3 Greedy algorithm 

We get an approximative solution to Problem A' if we set x* = in 
(fTTjl . This technique is called the greedy algorithm in Martello and 
Toth [8]. The resulting solution clearly satisfies the cost constraint 
due to the definition of the critical element. Its solution value 

s-l 
i=l 

satisfies 

In other words, when using the greedy algorithm, the error is bounded 
above by ps. 

We can expect the greedy algorithm to work well in financial ap- 
plications, since the probability space is usually relatively large, i.e. 
the probability for a single uJi is small. For example, in a binomial 
model with steps, the maximal probability for an individual u! is 
Pmax = max{p, 1 — p}^ , where p is the probability for an upward move 
during a single period. Even if p would be as high as 0.9, we would 
still have that, for instance for N = 100, Pmax would be as small as 
« 2.7 X 10"^. 
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Note that we can in general reach higher probabihty by truncating 
the critical Xs in (fTTI) than if we would truncate the entire set L := 
{dF/dQ = c* ■ E*{H)} in jT]). Indeed, after choosing all the elements 
in {dF/dQ > c* ■E,*{H)} we could still afford to choose some elements 
in L. 

However, from the computational point of view, it is often favor- 
able to group together the elements having the same pj/(?j-ratio. For 
example, consider a binomial model of A'^ steps. Then consists of 2^ 
elements. If H is an option whose value depends only on the value of 
the underlying asset at maturity, the ratio Pi/qi has the same value for 
all paths that lead to the same asset price at maturity. Thus, instead 
of ordering the 2''^ elements separately, we may order the pt/qi-leveh, 
whose number is at most + and search the critical level L. Finally, 
we can study the elements in L separately to see how many of them we 
can afford to hedge against. As a result, we obtain the same solution 
as with the greedy algorithm, but with less computational effort, since 



4 Maximizing the success probability when admissibil- 
ity is not required (Problem B) 

Next theorem shows that we can with any initial capital vq < K*{H) 
always create a so called quasi-replicating strategy, in other words a 
self-financing strategy that replicates H for all except one uj' . This 
fact is then used to give a solution to Problem B. The notion quasi- 
replicating strategy is discussed in context of binomial model e.g. in 
Favero and Favero and Vargiolu [3]. 

Theorem 4.1. Take vq < K*{H) and an arbitrary uj' G $7. Then we 
can construct a strategy ^ whose value process V satisfies Vq = vq and 
Vxioj) = H{ijj) for all oj ^ VL \ uj' . Moreover, 

vo-W{Hl^\^,) 

yT{uj) = — — — — . 



Proof. Note that Vt{uj') can be negative. However, due to complete- 
ness, there is a self-financing strategy ^ with value process V that 
replicates the random variable 

H - i^lov + 
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(See Remark EH). Moreover, Vq = E*(Vt) = W{H') = vq, since V is 
a P*-martingale. □ 

Corollary 4.2. Let vq <v and uji such that 

P(wi) = minP(w). 

Q 

Then the replicating strategy ^ for the random variable 
where 

P*(w*) 

is a solution to Problem B. 



5 Minimizing expected shortfall under admissibility 
condition (Problem C) 

The following result is a simplified version of Theorem 8.10 in Follmer 
and Schied [6], p. 341, where the result is proved for a general loss 
function and in a case where the market does not have to be complete. 

Theorem 5.1. Assume that there is a randomized test ip^ € IZ that 
maximizes the expectation K{Hip) among all ip G TZ satisfying the con- 
straint W{H ■ ip) < v. Then the replicating strategy ^ of the option 
H'-^ := H ■ ip^ solves Problem C. 

An optimal randomized test can be found by using the generalized 
Neyman-Pearson lemma, as is done in Follmer and Schied [6], p. 347. 
An alternative approach is to consider the problem of finding V'*^ as a 
continuous knapsack problem. Recall dlj and ^ and define 

™,:=MM:=«-i. ,12) 

Then the problem of finding an optimal randomized test can be written 
in the following form. 



Problem C. Find an n-dimensional vector that 

n 
i=l 



maximizes 
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under constraints 



Y^QiXi <a:= ^^^jjy < Xi < 1, i = 1, . . . ,n. 

This time the array {uji, . . . ,ujn) is ordered so that the quotient 
rrii/qi is non-increasing. Equations Q and (fT2]l give that 

rrii E*{H)pi 



where pi = P(a;j) and p* = P*(a;j). Thus, the states are ordered so that 
the quotient Pi/p* is non-increasing. The critical element t is defined 
as 



t := min \j-^qi> 



a 

I i=i ) 

and an optimal solution is by Theorem 2.1 in Martello and Toth [8] 
p. 16 given by 



X? 



1, i = l,...,t-l 

i = t 

^ 0, i = t + 1, . . . ,n. 



Note that does not in general coincide with x* in (fTT]l since the 
array (wi, . . . ,u;„) is ordered differently. 



6 Minimizing expected shortfall when admissibility is 
not required (Problem D) 

Favero [2] and Favero and Vargiolu [3] study the problem of minimizing 
expected shortfall in the special case of a binomial model, when admis- 
sibility is not required. Scagnellato and Vargiolu [H] discuss the same 
problem in a more general multinomial model. In those papers, the 
authors prove their results via dynamic programming. In this paper 
we provide an alternative approach by showing that even this problem 
can be reduced to a knapsack problem. To begin with, we state the 
following theorem. 

Theorem 6.1. Assume that the random variable maximizes E(X) 
among all random variables X that satisfy X < H and E*(X) < v. 
Then the replicating strategy ^ for X^ solves Problem D. 
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Proof. Take any self-financing strategy ^ witli value process V such 
that Vq < V and Vt < H . The random variable Vt satisfies E*(Vr) = 
Vo < f by the martingale property. Thus, we have by assumption that 

¥.{Vt) < E(X^). 

On the other hand, for the strategy .^^ with the value process we 
have 

and 

= ^*{V^) = E*(X^) < V. 

Further, the maximal expectation is attained by using strategy 
since 

W.{V^) =E(X^). 

□ 

Remark 6.2. Note that we cannot use in Theorem \6.1\ a similar ap- 
proach via randomized tests that is used in Theorem \5.1{ In Problem C 
the optimal solution is searched among strategies whose value processes 
V satisfy < Vt < H , which enables us to express Vt as the product 
of the claim H and a randomized test ■0- In Problem D, however, the 
value process Vt may become negative. On the other hand, the more 
general approach in Theorem \ 6.1\ could be used to reduce Problem C. 
Instead of searching for an optimal randomized test we could look for 
an optimal random variable X'~' that maximizes E(X) among all ran- 
dom variables X satisfying the constraints < X < H and¥,*{X) < v. 
In other words, the only difference would be the additional constraint 
X > 0, which is connected with the admissibility condition. 

Denote pi := p* := IP*(wj) and hi := H{uji). The problem of 

finding X^ obviously takes the following form in our finite probability 
space: 

Problem D'. Find a vector that maximizes 

n 

y^^xiPi 
1=1 

under constraints 

n 

XiPi < V, Xi<hi, i = l,...,n. 

1=1 
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This is almost similar to the continuous knapsack problem C. The 
only difference is that instead of constraint < Xj < 1 we have Xi < hi, 
i.e. the values of the decision variables Eire unbounded below and 
bounded by a deterministic, but varying non-negative boundary above. 
The problem can be solved in a way that resembles the solution method 
for a continuous knapsack problem. 

Assume that the array of the states Wj, i = 1, . . . , n is ordered so 
that 

Ply yPn 
P*l~"' ~ P*n 

The following theorem gives a solution to Problem D'. 
Theorem 6.3. An optimal solution x^ to Problem D' is given by 
xf = hi for i=l,...,n — 1 

Pn 

Proof. We prove our result in a way that resembles the proof of The- 
orem 2.1 in Martello and Toth [8]. Firstly note that for any optimal 
solution X it has to hold that 

n 

Y^Xip*=v. (14) 

i=l 

Without any loss of generality we can assume that Pi/p* > Pi+i/p*^]^ 
for all i. Let x* be an optimal solution to Problem D' and suppose 
that < hk for some k < n. Now if we take e > small enough, 
we could increase x^ by e and decrease x* by epl/p"^. But this would 
increase the value of our objective function by e{pk — PnPl/Pn) (> 
since Pk/p*k > Pn/Pn) give us a contradiction. Therefore, x\ = hk 
ior k < n is necessary for an optimal solution x* . The statement (fT3]l 
follows from (Ull) . □ 



Thus, we have proved that if we choose LOi such that 

dF , , dF , , 

— — (Wj) = mm — — u; , 
dF*^ ' uj&n dF*^ ' 

then the optimal strategy solving Problem D is by Theorem 16.11 the 
replicating strategy ^ for the random variable 
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where 

Note that this result coincides with the result in Scagnellato and Var- 
giolu [H], pp. 148, 149, where the authors prove it for a complete multi- 
nomial model via dynamic programming. 

The knapsack problem approach helps us to see the difference be- 
tween the optimal strategies when admissibility is or is not required, 
i.e. the difference between the solutions to Problems C and D. In both 
cases we arrange the states similarly, such that the ratio P/P* is non- 
increasing. Then we choose consecutively lo's for which we want to 
hedge the option perfectly, starting from uji that has the largest ratio. 
When admissibility is required, the critical element ujt is the first one 
that we no longer can afford to choose. The remaining capital is then 
used to construct a partial hedge for this state. If the admissibility 
condition is dropped, we can continue choosing w's until we reach the 
last one, LOn, having the least ratio. The value H^{ujn) is then adjusted 
so that the cost constraint is satisfied. 
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